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Abstract 

The distinctive features of passing the two-component extremely short pulses through the nonlinear media 
are discussed. The equations considered describe the propagation in the two-level anisotropic medium of 
the electromagnetic pulses consisting of ordinary and extraordinary components and an evolution of the 
transverse-longitudinal acoustic pulses in a crystal containing the paramagnetic impurities with effective spin 
S = 1/2. The solutions decreasing exponentially and algebraically are studied. 
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1 Introduction 

In recent years, an important subject of the theoretical studies became the coherent phenomena in the anisotropic 
media. This is caused by the significant progress in the technologies of growing the semiconductor crystals and 
producing the low-dimensional quantum structures. Unlike to the case of isotropic media, the parity of the 
stationary states of the quantum particles located in the anisotropic medium is not well defined. As a result, the 
diagonal elements of the matrix of the dipole moment operator and their difference usually called the permanent 
dipole moment (PDM) of the transition may not vanish. 

The dynamics of one-component extremely short electromagnetic pulses in the presence of PDM was consid- 
ered in works 1-5 . The stationary pulse solutions decreasing exponentially and algebraically of the full Maxwell- 
Bloch equations for the two-level medium were found 4 . It was shown that these solutions possess an asymmetry 
on the pulse polarity. The complete integrability of the corresponding system of the reduced Maxwell-Bloch 
equations with the help of the inverse scattering transformation method 6-8 was revealed in paper 2 . The multi- 
soliton solutions were constructed on fixed background, which was selected in such a manner that the problem 
under consideration is reduced to isotropic case. The effects of a pump on dynamics of these pulses were dis- 
cussed in 3 . Also, the numerical study of the pulse formation in the medium possessing PDM was performed 5 . 
An existence of solitary stable bipolar signal with nonzero time area was discovered. The asymmetry on the 
polarity manifested itself there in that the sign of the time area is determined by the sign of PDM. 

The propagation through optically uniaxial medium of the two-component electromagnetic pulses consist- 
ing of short- wavelength ordinary and long- wavelength extraordinary components was investigated in 9-12 . As 
opposite to the case of one-component pulses mentioned above, an application of the slowly varying envelope 
approximation does not eliminate the influence of PDM in this problem. It was shown that, under conditions 
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of a strong interaction between the components, the regimes of the pulse propagation may differ from the self- 
induced transparency. Various regimes of the optical transparency were classified with respect to the pulse 
velocity and the degree of the resonant medium excitation. The two-component pulse solutions possess also the 
PDM-induced asymmetry: signs of the pulse extraordinary component and PDM are opposite. It is remark- 
able that the equations distinguishing from ones studied in 9-12 by notations only describe the two-component 
(transverse-longitudinal) acoustic pulse propagation in the low-temperature paramagnetic crystal 13,14 . 

To investigate the propagation of the two-component pulses in optically uniaxial medium, we apply here 
the approximation of the unidirectional propagation instead of the slowly varying envelope approximation. The 
two-component system of reduced Maxwell-Bloch equations is derived and, as in the case treated in works 9-12 , 
its acoustic analogue is found. The properties of exponentially and algebraically decreasing pulse solutions of 
these equations are discussed. 



2 Two- component system of the Maxwell— B loch equations for an 
anisotropic medium and its acoustic analogue 

Consider the optically uniaxial medium containing the two-level quantum particles. Let the medium optical 
axes coincide with z axes of the Cartesian coordinate system. Suppose that the plane electromagnetic pulse 
propagates through the medium in the positive direction of y axes. Then, ordinary E Q and extraordinary E e 
components of the electric field of the pulse are parallel to x and z axes, respectively. The Maxwell equations 
yield in the case we study the next system for the electric field components: 
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where P Q and P e are the ordinary and extraordinary components of the polarization, which is induced by the 
two-level particles; n Q and n e arc ordinary and extraordinary refractive indices, respectively; c is the speed of 
light in free space. 

To describe the evolution of the state of the quantum particle, we make use the von Neumann equation on 
density matrix p: 

ih^ = [H,p]. (3) 

Here Hamiltonian H of the two-level particle is defined as given 

H = diag(0, hu ) - d x E a - d z E e , (4) 

where u>o is the resonant frequency of quantum transition 1 1)— >|2); h is the Plank's constant; d x and d z are the 
matrices of the projection of the dipole moment operator on x and z axes, respectively. 

In the geometry chosen, the expressions for the polarization components are 

P = NTr(pd x ), (5) 

P e =NTr(pd z ), (6) 

where N is the concentration of the quantum particles. Since the medium is axially symmetric, the matrices of 
the dipole moment have the following form: 
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Then Hamiltonian Q reads as 

»-( 

-di 2 E hui - D22E, 
and equation J3J) gives the next system for the density matrix elements: 
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where D — D\\ — D22 is PDM of the transition. Formulas 0), (0) and ©-© imply 
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Let us suppose that concentration of the quantum particles is small (irdi 2 N/(hujo) <C 1)- The order of 
derivatives in the equations for the electric field components can be reduced in this case. Indeed, by applying 
the approximation of the unidirectional propagation 15 to equations 1|12(1 . (|13(l and excluding the derivatives of 
the elements of the density matrix with the help of formulas illOt and JUJ, we obtain: 
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It is convenient for the subsequent consideration to rewrite equations (|10fl , and (|14|) , l|15|) in dimensionless 
form 
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Obviously, system l|l(j|) - (|19[l coincides with the reduced Maxwell-Bloch equations for isotropic medium if 
k = v = 0. One can see that the electric field components fulfill different functions here. The ordinary component 
causes the quantum transitions, while the extraordinary one shifts its frequency. Also, the same functions were 
executed by these components in works 9-12 , where the system, which follows equations ljl 6|) - (|1 9|) in the slowly 
varying envelope approximation, was considered. 

ft is remarkable that the system of equations equivalent to i|l 6(l - (|f 9|l appeared in the acoustics. Namely, 
the propagation of the transverse-longitudinal acoustic pulses in a crystal containing paramagnetic impurities 
with effective spin S — 1/2 in a direction parallel to the external magnetic field is described by next equations 
(see 16 - 17 ): 
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E and U are dimensionless variables defining transverse and longitudinal components of the acoustic pulse, 
respectively; S x , S y and S z are expressed through the elements of the density matrix of the paramagnetic 
impurities; x an d V are dimensionless spatial and retarded time variables; constant b is defined through the 
coupling constants of the spin-phonon interaction; v\ and V2 are linear velocities of the transverse and longitudinal 
acoustic waves, respectively; u>b is the frequency of the Zeeman splitting of the Kramers's doublets; n is the 
concentration of the paramagnetic impurities; no is the mean crystal density; g is the component of the Lande 
tensor; / 3 is the coupling constant of the spin-phonon interaction. 

The variables of both systems are connected by relations: 

E = ku, U=J— (1 + v), (26) 

V n o 

S x = -k(a + a*), S y = ik(a-a*), S z = 2ka 3 , (27) 

ii=jt. (28) 

(29) 

Complete integrability of equations l|2i p — l|25 |l in the frameworks of the inverse scattering transformation 
method 6-8 in the case, when the linear velocities of transverse and longitudinal acoustic waves are equal, was 
established in paper 17 . Consequently, system lfTH|l - (|Tl^l is also integrable if 6 n = 0. 
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3 Simplest solutions 



Let us begin with system I|16|) - (ll9fl . At first, we suppose that the refractive indices of the medium are equal: 
n e — n a . As it was noted at the end of previous section, system we deal with is integrable by inverse scattering 
transformation method in this case, and its multisoliton solutions can be obtained by means of the algebraic 
methods of the soliton theory. The Darboux transformation 18 technique can be used for constructing these 
solutions, for instance. Since the one-soliton solution is stationary, it can be found by direct integration of 
equations (|16[) ~H19 [) . After straightforward calculations, one reveals that there can exist two families of the one- 
soliton solutions. These families arc distinguished by the domains on the spectral parameter plane, in which the 
points of the discrete spectrum of the solution spectral data lie. Solution of the first family has next form: 
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(Jo is an initial population of the quantum level (|<to| < 1/2). Free parameter of this solution is real constant T. 
It is assumed here that parameter T takes the values such that condition A > holds. This condition is true 
with arbitrary T, if K > 1. For K < 1, parameter T has to satisfy constraint 

9 K 
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In the case A < 0, the one-soliton solution is written as follows 
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This solution belongs to the second family of the one-soliton solutions. It can exist only if K < 1. Also, the 
parameter T value is such that 
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The well-known one-soliton solution of the reduced Maxwell-Bloch equations 15 is obtained from formulas Q34JI - 
(j23) if we put k = 0. 
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It is seen from equations <|2U|) , (|31|) and (|35|) that there exists asymmetry on the polarity of the extraordinary 
component of the pulse electric field: signs of E e and PDM are opposite. Similar asymmetry was found in 9-12 
for two-component pulses, in which the ordinary component has slowly varying envelope. The number of peaks 
of variable u of the second family solution depends on parameter T. If 



T 2 < 



2K 



l-K 

it has two peaks, while at 

2 

T > 



1 -K 



only single peak exists. One can show that a degree of the medium excitation grows with increase of \T\ for the 
first one-soliton solution family and, on the contrary, with decrease of |T| for second one. The full inversion of 
the quantum level population can take place only if K < 1. The curves of variables u and cr 3 of the first and 
second families of one-soliton solution are presented on Figs 1, 2 and Fig. 3, respectively. The distance between 
the peaks of u depends on parameter T and can take an arbitrary value if K < 1 (see curves on Figs 2a and 3a) . 
If K > 1, there exists an upper limit of this distance that corresponds to Fig. 6 below (compare with Fig. la). 
The strongest excitation of the quantum particles occurs in the center of the pulses in all cases. 




Fig. 1. Profiles of u and 0-3 of the first one-soliton family with K = 1.05, 00 = —0.5, T — 2 (solid curves) 
and T = 10 (dotted lines). 



Consider different limiting cases of the one-soliton solutions. If K < 1, we can tend A to zero. Then, after 
appropriate shift alone r (or £) axes, we obtain from formulas (|3"0|) (|33|l or (|3^|) - (|3Tjl following expressions: 
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It is necessary to put T — k/\/A — k 2 in this solution in accordance with the limiting procedure. If |T| — > 00, 
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Fig. 2. Profiles of u and 03 of the first one-soliton family with K = 0.8, <ro = —0.5, T = 1.4 (solid curves) 
and T = 1.9999 (dotted curves). 




then formulas l|3Ufl - 133 (1 lead to algebraic solution: 
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This limit is allowed in the case K > 1 only. 

Dependence of variables u and 03 of solutions (|38[) - l|41[) and (!4Lj|> — (|45|> on r are presented in Figs 4, 5 and 6, 
respectively. Comparing with curves on Figs 2a and 3a, we can say for the variable u case that the first limiting 
procedure (A — > 0) corresponds to moving one of its peaks into infinity. The pulse, whose velocity exceeds the 
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Fig. 5. Profiles of u (solid curve), v (thick curve) and 03 (dotted line) 
in the case A = with K = 0.8, T = -2.0, cr = 0.5. 



light velocity, is plotted in Fig. 5. Note that variable v tends to —2 in the limit t — > —00. For this reason, the 
energy levels trade their places owing to the Stark effect. 

Let us now consider the case of unequal refractive indices (n e 7^ n e ). If we define coefficient K as 

n k 2 ( 1 + T 2 



in the formulas (JSOJ-ljSSl) and lp^jl - (P7jl . then they give solutions of system l[Tfi|l - (fB^I with S n ^ 0. Since 
coefficient K characterizing the anisotropy of the medium depends now on the parameter of the pulse, we can 
say that the anisotropy becomes effective. It is seen that this coefficient is unbounded if 

rp2 _ ' 



1 + 4a /5 n 



and it can change the sign. Similar results were obtained in works 9 12 , where the approximation of the slowly 
varying envelope was applied for the description of the dynamics of the two-component pulses. 

Obviously, the limiting cases studied above can be also modified to satisfy equations (|lfi|) - H19|) with n e 7^ n a . 
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Fig. 6. Profiles of u (solid curve) and 03 (dotted line) of the algebraic 
solution with K — 1.15, ao — —0.5. 



In particular, one has to put 



n a k 2 ( S n 
K = — —n e I 1 + - — 

4(Tn 



in formulas (|42|I - H45|I . Finally, coefficient K should be determined by relation (|46() in equations (|38[) - i|41[) . while 
parameter T should be equal to 



where 



se — — k — 
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The solutions of system <|21|) (I25fl can be found using its equivalence to system (|16|1 - <|19[) (see relations 
1261) (|29Jl l. In particular, formulas H30(l - (|33[1 and H34(l - (|37|l give two families of the pulse solutions of the acoustic 
system: 
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So is an initial population of the quantum level, real constant T is free parameter of these solutions. Coefficient 
B has to be positive for the former solution and negative for latter one. For the sake of simplicity, we suppose 
here that S v = 0. The main properties of the pulse solutions presented are, obviously, the same as for the system 
H16|) (|19|) solutions. The appropriate modification of U and b in these formulas allows one to obtain the solutions 
of system (JHJ-lESJl with S v ^ 0. 



4 Conclusion 



The two-component system of the reduced Maxwell-Bloch equations for anisotropic medium and its acoustic 
analogue have been considered in this report. As it was done in works 11 ' 12 , the linear velocities of both the 
components are not supposed to be equal. Two families of the one-parameter extremely short pulse solutions 
and their limiting cases, including one that leads to algebraically decreasing two-component pulse, are studied. 
The ordinary component of the pulse solutions found can have one peak like in the isotropic case or, also, two 
peaks, whose signs can be identical or opposite. If PDM is small (\D/2d\ < 1), then the distance between the 
peaks can be unbounded. Nevertheless, the strongest degree of the medium excitation is always achieved in the 
center of the pulses. The asymmetry on the polarity of the signal takes place here similarly to that in the cases 
of one-component pulses ' 5 and two-component pulses, in which the ordinary component has higher-frequency 
filling 9-12 . Namely, the signs of the extraordinary component and PDM are opposite. 
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